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ABSTRACT 


A  variant  of  the  simplex  method  is  given  for  solving  linear 
programs  with  M  +  L  equations,  L  of  which  have  the  property  that 
each  variable  has  at  most  one  nonzero  coefficient.  Special  cases  Include 
transportation  problems,  programs  with  upper  bounded  variables, 
assignment  and  weighted  distribution  problems.  The  algorithm  described 
uses  a  working  basis  ^f  M  rows  for  pi  voting,  pricing,  and  inversion 
which  for  large  L  can  result  in  a  substantial  reduction  of  computation. 
This  working  basis  is  only  M  x  M  and  is  a  further  reduction  of  the  size 
found  in  an  earlier  version,  see  [lj.  Unfortunately,  to  achieve  this 
reduction,  row  as  well  as  column  transformations  must  now  be  made. 


GENERALIZED  UPPER  BOUNDING  TECHNIQUES  -  II 


by  G.  B.  Dantzig  and  R.  M.  Van  Slyke 


1.  Introduction 

A s  in  [l]  we  assume  that  we  are  concerned  with  a  linear  program 
in  which  each  variable  has  at  meat  one  non- zero  coefficient  in  which 
is  the  last  L  equations  to  be  nonnegative  and  the  corresponding  last 
L  constant  terms  are  positive.  In  section  IV  we  indicate  the  necessary 
modifications  to  handle  negative  coefficients  in  those  equations.  By 
normalizing  the  variables  and  multiplying  the  equations  by  constants,  we 
can  assume  without  loss  of  generality  that  all  nonzero  coefficients  and 
constants  in  the  last  L  equations  are  ones,  (l). 
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Ve  review  here  for  completeness  some  necessary  definitions  and  theorems 
we  will  borrow  from  [1], 

The  L  set  of  variables  or  columns.  ,  will  refer  (depending 

on  context)  to  those  variables  or  columns  corresponding  to  the  columns 

of  coefficients  In  (l)  with  1  as  their  M  +  component.  SQ  ,  the 
th 

0  set,  is  the  set  corresponding  to  columns  with  zeros  for  the 

st  th 

M  +  1  through  M  +  L  coefficients. 

We  assume.  that  the  system  (l)  is  of  full  rank  and  denote  by 


^1  ^1  0 
[A  ,...,A  ]  ,  a  basis  for  the  system.  We  always  assume  A  ■  A 

the  coefficient  of  the  variable  to  be  optimized.  The  underscoring  is 

to  differentiate  coefficient  vectors  with  all  M  +  L  components  from 

the  reduced  vectors  of  the  first  M  coefficients  which  -ill  be 

denoted  without  the  underscoring.  There  will  be  no  underscoring  for 

individual  components  of  these  two  different  types  of  vectors 

even  though  they  differ  in  the  number  of  their  components. 


Theorem  l)  At  least  one  variable  from  each  set  (with  the  possible 

exception  of  SQ  )  is  basic. 

Theorem  2)  The  number  of  sets  containing  two  or  more  basic  variables  is 
at  most  M  . 


The  sets  containing  two  or  more  basic  variables  plus  SQ  are  called 
essential  sets  (the  term  here  is  used  in  a  slightly  different  sense 
than  in  [1]).  An  essential  set  for  one  basis  may  become  an  unessential 
one  in  the  next. 

In  the  next  section  we  outline  the  method,  in  the  following  we 
formalize  it  as  an  algorithm.  In  the  last  section  we  indicate  the  modifi¬ 
cations  required  to  handle  negative  coefficients  in  the  laat  L  equations 
and  finally  in  the  appendix  the  method  is  carried  out  on  an  example. 


II.  The  Method 


Given  a  feasible  basis  ,  ve  assume  we  have  selected  for  each  S. 

k  1 

cue  basic  variable  x,  to  be  the  key  variable.  A  is  said  to  be  the 
key  column.  For  the  key  column  can  be  taken  for  convenience  to 
be  a  dummy  column  with  all  zero  components.  We  then  consider  the 
system  obtained  by  subtracting  the  key  columns  from  every  other  column 
in  their  respective  sets  (in  (2)  ve  assume  for  simplicity  that  the  key 
variable  was  the  first  one  in  each  set).  In  this  modified  system  the 
value  of  the  key  variables  must  clearly  be  one.  We  treat  these  variables 
as  ve  vould  variables  at  upper  bound  in  an  upper  bounded  variables 
algorithm  for  the  simplex  method  and  subtract  their  coefficients  from  the 
right  hand  side. 


D  1  rn  A  1 


J 


AJ  - 


J  i 


(3) 


d  .  b  - 


i-0 


1  Obtaining  a  first  feasible  solution  is  accomplished  using  this  method 
vlth  a  phase  I  set  up  as  in  the  usual  simplex  method. 
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where  the  last  L  components  of  d  are  zero. 

We  assume  that  for  j  »  (key)  to  be  absent  from  the 

system.  The  working  basis.  B  ,  is  given  by  B  »  (D^  |  is  basic  and 
not  key) .  By  virtue  of  theorem  1  it  is  clear  that  B  has  exactly 
M  columns.  We  assume  B1  ■  corresponds  to  the  coefficient  of 
the  variable  to  be  optimized.  We  define  the  derived  system  to  be 


5>/ 


*  d  and  it  is  easy  to  prove: 


Theorem 


B  is  a  basis  for  (4). 


Proof :  Suppose  IV-0-  Since  B^  differs  from  B“  (the 

same  column  considered  in  the  system  depicted  in  (2))  by  only  0 
components  ^  ®  •  But  this  implies  that  the  plus  the 

key  columns  are  linearly  dependent  since  the  by  themselves  are 

linearly  dependent.  On  the  other  hand  this  set  is  obtained  from  a 
(nonsingular)  basis  by  subtraction  of  columns  from  within  the  set 
which  does  not  reduce  the  rank,  yielding  a  contradiction. 

By  Theorem  2  there  exist  at  most  M  sets  with  more  than  one 
basic  variable.  These  sets  and  SQ  are  the  only  sets  which  contain 
members  of  B  ;  i.e.,  they  are  the  essential  (including  SQ)  sets. 

Thus,  with  each  feasible  basis  for  the  original  system  (l),  we 
have  associated  a  set  of  L  key  variables  (and  one  dummy)  and  a 
basis  for  the  derived  system.  We  now  show  that  we  can  carry  out  the 
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steps  of  the  simplex  method  using  Just  the  inverse  B  1  of  B  ,  the 
reduced  basis,  and  the  corresponding  basic  solution,  d  ,  of  the 
derived  system  (4). 

Hie  first  step  is  to  obtain  a  set  of  prices  for  (l).  Let 
us  denote  by  77  -  (7^,. . .  ,7r^)  the  prices  on  the  first  M  equations 
and  |i  ■  the  prices  on  the  last  L  .  These  prices  are 

determined  uniquely  by  the  condition  that 

J, 


(tt  ,  p.)A  -  (tt,  *i)A  1  -  1 
(tt  ,  n)A  -  0 


i«2,. . . ,M+L 


A  m"\ 

Let  7T  ■  (B  )  ,  the  first  row  of  the  inverse  of  the  working  basis 

B  .  It  has  the  property  that 


ftB1  -  ttA°  -  1 


ttB' 


J 


J*2,. . .  ,M 


i.e. ,  tt  is  a  set  of  prices  for  (4).  To  extend  it  to  a  set  of 
prices  for  (2)  is  trivial  we  simply  set 


(5) 


A  *  A  JL 

M-s"  -  7tA 


/®1 ,  •  •  • ,  L  t 

J, 


Now  consider  for  basic  columns  A 


J,  h  J 

(tt  ,  (i)A  1  -  (  tt,  Cl)A  1  m  o  if  A  1  is  key  (5) 


or 


4 

%j  j 


(tt  ,  ft)A  1  ■  (tt  ,  £)(§?■  A*)  for  some  k 

-  fir  ,  (Db1  ♦  (tt  ,  (l)Ak 
^i 

■  0+0  if  A  la  not  key  . 


Thus  (tt  ,  M- )  is  a  set  of  prices  for  the  original  system  (l). 
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T 


't' 


Using  these  prices  we  can  "price  out"  the  columns  of  (l)  to  find 

one  to  enter  into  the  basis.  Using  the  usual  simplex  criterion,  the 

s 

incoming  column  A  would  be  chosen  by 


where 


(if  i  | jl)A8  ■  min  (tt  ,  p.)A^ 

J 


min  A 

J 


i 


Vi +  ^ 


for  A°  €  S, 


If  A  >  0  ,  we  have  an  optimal  basic  feasible  solution  and  we're 

s  ® 

8  8 

done;  otherwise  we  bring  A  into  the  basis.  Assume  A  €  3^  . 

s 

To  do  this,  we  must  express  A  and  b  in  terms  of  the  current 
basis  for  (l).  If  we  let 

D  8  .  B-1DS  .  B_1(AS  -  A  °) 


then 

M 

(6)  (A*  -  A S  .  “b1  .  ][d  “(A^  -  AVl) 

i.1 

where  indicates  the  column  number  in  (2)  corresponding  to  the 

i  column  of  the  working  basis  and  denotes  the  column  number 

of  the  corresponding  key  variable. 

s 

We  denote  the  representation  of  A  in  terms  of  the  current 
basis  by  A  ^  ,  that  is: 


\ 

U  j 


M+L 

I*;*'1 

i-1 


Fran  (6)  we  see 
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(7) 


A  8 


1  - 1 
Vt-kB 


if  A 


D  *  if  A1 


•I* 

V-l 


~  V  if  A± 


nt 

A  for  some  t 


A  for  some  t 


0  otherwise 


The  current  values  for  the  variables  in  the  basis  b^  are  given 
either  by  updating  the  values  of  the  previous  iteration  in  the  usual  way 

or  recomputed  in  a  similar  way  to  A  above.  That  is 

J 


(6) 

then 


Let  d  -  dj^. . .  ,dM  be  given  by 
d  «  B_1(b  -  ^  A  l)  .  B_1d 


(b 


and  as  in  (7)  the  b^  are  given  by 


(9) 


1  - 


if  A  x  Is  key 


bi  ■< 


VJi  J  TJt 

d  if  A  ■  A  for  some  t 

w 


0  otherwise 


Finding  the  variable  to  leave  the  basis  is  accomplished  in 
exactly  the  same  way  as  in  the  ordinary  simplex  method.  Let 
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(10) 


b 


■  min 
>  0 


b 

i 


,  i  ®  y  .  .  .  y  M  ^  L  . 


_  S  j 

where  we  of  course  require  that  A  >  0  .  Let  us  assume  A  r  e  S  . 

r  P 

Three  cases  can  occur  in  the  updating  process : 

L 

(a)  If  is  not  essential  and  A  6  Sfl  ;  i.e. ,  the  outgoing 

variable  is  the  key  variable  in  S  then  B  remains  unchanged, 

j  0 

and  A8  simply  replaces  A  r  as  the  key  variable  in  . 

This  requires  the  updating  of  d  which  is  accomplished  as 
follows : 


-  d  -  b’V  .  kT) 


d  -  D  8 


J  k 

Observing  that  A  r  =,  a  0  ,  this  is  easy  to  compute  since  we 
already  have  d  and  the  second  term  was  generated  in  determining  the 

i;. 

^r  -1 

(b)  If  A  is  not  a  key  variable,  then  we  update  B  simply  by 


—  s 


J. 


pivoting  on  the  column  D  on  the  row  which  A  -  A  M  occupies 


1 J  Where  the  symbol  *  "  does  not  indicate  equality  but  rather 

that  the  expression  on  the  right  replaces  or  (updates)  the 
variable  on  the  left. 


-9- 


in  the  working  basis.  In  symbols  B  ^  :  *  PB  1  where  P  is 
the  matrix  which  performs  the  pivot,  d  is  updated  by  applying 
this  pivot  to  the  old  d  . 

J 

(c)  If  A  r  €  S  is  a  key  variable  in  a  different  set  than  the  one 
P 

8 

containing  A  €  ,  we  must  first  change  the  key  variable  in 

S  .  To  change  the  inverse  of  the  pseudo  basis  B  ,  we  consider  all 
p  i 

J  Jr 

columns  of  B  of  the  form  A  -  A  there  must  be  at  least  one 

such  since  after  A8  enters  the  basis  S  must  contain  a 

P 

basic  variable  (Theorem  l).  One  of  these  columns,  call  it  A  , 
is  to  become  the  key  variable.  To  get  the  new  pseudo-basis  B 

from  the  old  one  B  we  wish  to  multiply  the  column  (in  B  ) 

k  Jj. 

A  -  A  by  -1  to  obtain  A  -  A  and  we  wish  to  subtract 

k  *^r  j 

A  -  A  from  every  other  column  of  the  form  A  -  A  for 

J  k 

J  j*  k  to  obtain  A  -  A  .  Hiat  is 


B  :  ■  BT  where  T  is  of  the  form 


(12) 


0...0-1...-1  -1  -1...  -1  0...0 

1 


Where  the  -1 '  s  occur  in  the  columns  corresponding  to  A'*  €  S  ,  and 
the  row  corresponds  to  the  new  key  variable,  A  . 


B 


-1 


-1-1 
T  B 


arid  it  is  easily  verified  that  T  =  T 
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T 


j* 

since  applying  the  process  twice  replaces  A  as  the  key  variable. 

The  values  for  d  are  updated  by  applying  T  ^  .  Now  with  the  new 

key  variable  in  S  we  simply  apply  the  process  outlined  in  b  . 

P 

With  our  updated  B  ,  y  and  key  variables  we  sire  now  ready  to 
make  another  iteration.  If  the  inverse  of  the  pseudo-basis  is  expressed 
in  product  form  we  have 


where  each  T*  is  either  of  the  form  (12 )  or  (13),  the  latter 

th  T 

resulting  from  a  pivot  on  the  r  element  of  the  column  (ol  , . . .  ,a  )  . 


(13)  P  - 


-V°r 


-a 


r-l 


-a  /a 

m  r 


1 


As  W.  Orchard-Hays  has  remarked  to  one  of  the  authors  we  can  if  we 
wish, express  each  transformation  of  the  form  (12)  as  a  sequence  of 
transformation  of  the  form  (13).  Suppose  we  wish  to  express  a  matrix  T 
in  terms  of  transformations  of  the  form  (13),  where  -l's  appear  in 
columns  hy  hy . . .  ,  h^  and  suppose  that  the  -1  in  column  h. 
lies  on  the  diagonal;  in  other  words,  all  the  -l's  are  in  the 

a 

h  u  row.  Let  P(r,s)  denote  the  pivot  matrix  (l 3 )  with  a,  ^  0 

U  J 
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J  f  r,  or  a  =1  and  a  ■  -1  .  When  multiplying  on  the  left 
r  s 

th 

this  matrix  has  the  effect  of  subtracting  the  r  row  from  the 

s  h  row.  Finally  let  P(r,r)  be  the  matrix  with  all  plus  ones  on 

th 

the  diagonal  except  in  the  r  diagonal  element  which  is  -1  . 

Every  other  element  is  zero.  When  multiplying  on  the  ri^ht . 

th 

this  has  the  effect  of  multiplying  the  r  column  by  -1  .  It 
is  then  easy  to  see  that 

T .  p(V  vp(>v  Vi>---p(vyp<v  v 
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FLOW  DIAGRAM  OF  THE  GENERALIZED  UPPER- BOUND  ALGORITHM 


(1) 


ENTER 


Enter  with  inverse  of  working  oasis 
B  \  uasic  solution,  d  .  .l  U-e 
derived  system  (U)  and  a  se  of  L 
kc\  ariaoles. 


U 


isl _  _ I _ 

|  Price  out  coLumns  usin^  (3  ^  ) 


bin  n  ^  O 


u 


Min 


=  u  <  O 


00 


J  ~S  €  S° 


—  a 


Find  D  (  a  ,  and  t  ,  and  use 

J  i' 

row  criteria  to  find  A  the 
variaole  to  oe  dropped.  _ 


►  - 


Jr 
H  £  S 


I- 


noo  xe„ 

- ■* 


Pivot  on  D  in  the  r 


ow 
r  . 


corresponding  to  D  in 
3  ^  and  update  d  . 


Return  to  (l)  for  next 
i  t.eration. 


.  P  t  0 

.Co) _ \ _ 

Mane  some  oasic  column 

in  set  S  otner  than 

-1 


A  r  ,,ev.  Update  B  ,| 
and  d  . 

. _ r. _ 


Go  to  (i?) 


lY-rminatc  cuirei. 


solution  is 

optimal. 


r 


(V) 

P  *  o 

Mu*  e  AJ  ,-pj  i  s  cu  i 
of  A'  1  and  upda  e 

h  . 

Return  to 

for  .L..tj 

i teration. 

1 

FIGURE  1 


III.  Description  of  Algorithm 


Referring  to  figure  1,  the  algorithm  takes  place  in  the  following 
steps : 

(1)  We  assume  we  enter  the  algorithm  with  the  inverse  B  ^  of  the 
working  basis,  the  value  d  of  the  appropriate  basic  solution 
of  the  derived  system  (4)  and  the  set  of  key  variables.  To  get 
this  initial  solution,  the  usual  phase  I  procedure  can  be 
carried  out  in  the  obvious  way. 

(2)  Let  7T^  -  (b"1)*  for  1  "  6111,1  for  each  8Ct  i  0) 


M 


let  p, £ 
Let 


V  £  / 

^Ai  *  where  A  1  is  "the  key  column  in  . 


i«l 


M 


Z  Vl  +  H 


for  A°  e  S 


/  * 


i-1 


Let  A  -  “in  A.  and  suppose  A  e  S  .  If  A  >  0  ,  we  go 
s  J  a  b  fc 

to  Step  (3);  otherwise  skip  to  (4). 


(3)  Terminate  we  have  an  optimal  solution. 


(4)  Find  D8  -  B’1 


-l/.s 


s 


D«B(A-A),A  by  means  of  Equation 
(7)  and  b  by  means  of  Equation  (9).  Use  the  usual  simplex 
decision  rule  Equation  (10)  to  find  the  variable  to  be  dropped 

Jr  Jr  Jr 

A  and  suppose  A  e  .  If  A  is  key,  go  to  step  (3); 
if  not  and  p  ^  0  ,  go  to  (6);  if  p  ■  0  ,  go  to  step  (7). 

/  \  ”8  r 

(3)  We  pivot  with  respect  to  D  in  the  rov  corresponding  to  D 
in  B  ^  and  update  d  by  applying  the  pivot  transform  to  it. 
We  then  return  to  step  (l)  for  another  iteration. 
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k. 

(6)  Make  some  basic  column  say  a  &  +  j  in  set  S  key 

J  -1  P  P 

instead  of  A  r  .  Update  B*1  by  applying  a  column  trans¬ 
formation  of  the  form  (12)  and  update  B  1  by  B  1  :  ■  T  J'B  \ 

d  is  updated  by  d  :  ■  T  ^d  .  We  then  can  go  to  step  (3). 
s  Jr 

(7)  Make  A  key  instecui  of  A  and  update  d  by 

d:«  d  -  D  8  . 

Return  to  step  (l). 

IV.  Negative  Coefficients : 

When  negative  coefficients  appear  in  the  last  M  +  L  equations, 
the  algorithm  is  changed  in  a  quite  obvious  way.  We  can  assume  without 
loss  of  generality  that  the  coefficients  in  the  last  L  equations  are 
+1  or  -1  and  the  last  L  right  hand  side  components  are  +1. 

Theorems  1,  2,  and  3  still  hold,  and  we  can  require  that  each  key 
column  have  a  +1_  in  the  last  L  equations  since  clearly  each  set 
must  have  such  a  column  which  is  basic.  In  the  pricing  process  if 
the  column  Ac  to  be  priced  has  a  negative  coefficient  in  the  last 
L  components,  the  appropriate  4  is  subtracted  rather  than  added 
to  7tA^  .  To  form  the  difference  columns  Dv  the  key  column  is 
added  to  columns  with  a  -  1  rather  than  subtracted  and  appropriate 
modifications  in  equations  (7)  and  (9)  must  be  made  to  reflect  this. 

Other  than  these  slight  modifications,  the  algorithm  proceeds  exactly 
as  before. 
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I  o'  I 


Consider  the  following  example  (Figure  2)  with  Ma  3 


Figure  2 

The  initial  basis  is  aV^A-Wa^  and  A^A^A®  are  key 
The  working  basis  B  =  {A°  -  0  ,  A?  -  ,  A  -  A"1"}  ;  hence 


which  has  an  inverse : 


With  the  aid  of  (5)  we  find  the  prices  [  7T ,  (j.  J  =  [4  >  4  >  4  > 

.2  _2  2] 

2  >  g  >  2  ^  c  ' 

We  then  price  out  and  find  A7  €  wins. 


A  J  -  A  6  «  b"V  -  A6] 


4  4  4” 

t  -t  -t 

001 


i.e.,  A7  -  A^  -  -3A°  -  4(A^  -  A^) 
or  A7  «  -5A°  +  ^AX  -  ♦  A^ 

7 

giving  a  representation  of  A  in  terms  of  the  full  basis. 
A7  -  (-3,  J,  -i,  0,  0,  i,  0)T 
We  obtain  the  values  of  the  variables  by  considering 

B_1  [b  -  r  A£)  -  B_1 

This  means 

b  -  A  1  =  5A°  -  JA^  ♦ 
or 

b  »  5A°  +  ^A1  +  4A2  +  A4  4  A  +  A  +  A 


hence 


We  now  determine  the  variable  going  out  of  the  basis  by 


9- min  ^  5 
Ilo  K  K 


and  r  could  be  2  or  7  .  Taking  it  to  be  7  ,  we  see  that  since  the 

6  7 

set  is  inessential  and  ©  ■  1  we  just  replace  A  by  A1  as 

a  key  variable  and  B  remains  unchanged.  The  new  multipliers  are 

r  Z  i  ri.  —  —  _Z.  2] 

L  •>  f  p.  J  *  ^2^22^  2  9  2  9  2  9 


and  this  time  A  €  prices  out  optimally. 


b'X(A3  -  A6] 


i  i  i' 

i  -i  -i 

0  0  1 


that  is,  A 


-A°  ♦  5A1 


2  8 

5A  +  A  and 


A*  -  [-1,  5,  -5,  0,  0,  0,  1,  0] 


-1  V  K£ 
B  lb  -  A1] 


i  i  i 
i  -i  -i 


-18- 


b  =  6A°  +  -  A1  +  A1  +  A4  +  A  +  A  +  A 

-  6A°  ♦  A2  +  A4  ♦  A>  ♦  A7  ♦  A8 

b  «  t  6,  0,  1,  0,  1,  1,  1,  1  ] 

h  ho. 

9  =  mj.n  m  r  »  r  *  2 

A/>0  Ai  A2  Jr  -  1 


therefore,  we  want  to  drop  f  which  however  is  key.  So  first 

1  2 

we  must  replace  A  by  A  as  a  key  variable.  To  do  this  we  take 
our  current  working  basis 

t,  r .  0  2  1  3  1, 

B  =  [A  -  0  ,  A  -  A  ,  A  -  A] 


and  postmultiply  it  by 


T  - 


10  0 
0  -1  -1 

0  0  1 


which  has  the  effect  of  subtracting  the  second  column  from  the  third 
and  reversing  the  sign  of  the  second  column. 


B' 


HT  = 


(b1  )  ~1  =.  t"1b"1 


where 


T 


-1 


10  0 
0  -1  -1 

0  0  1 


Hence 


-19 


(S') 


-1 


1 

0 

0 


0  0 

-1  -1 


i  i  i 

i  -i  -i 


i 

l 

T 


i  i 

i  4 


We  then  pivot  in  the  vector  column 

i-l  1 

2  2 


/  ,x-l/  9  8X 

(B* )  (A  -  A  ) 


1 
4 

0  0 


1 

T 


1 

2 

J. 

4 


'-11 

9 


-1 

5 

0 


on  the  second  component. 

This  gives  us  a  new  inverse  basis 


B 


-1 


PB' 


-1 


rl  i. 

5 

0  J 


0 

0 

1 


9  11 
20  So  20 

i.  i_-i_ 

20  20  20 


1 

2 

1 

'4 


1 

2 

1 

4 


1 

2 

4 


20- 


T 


The  new  prices  are 


,1  i  L.l  1  1  «) 

'20  ’  20  ’  20  '  20  ’  20’  20  ’  20  ’ 


and.  upon  pricing  out  we  find  that  all  columns  price  out  nonnegatively 
and  the  optimal  solution  is  given  by 


B 


2 

11 

L_ 

I 

— i  k 

20 

20 

20 

1  [b  -  1 

1 

1 

L. 

20 

20 

‘20 

0 

0 

1 

6 

6 


6 

0 


L°  j 


and 


2  4  5  7  8  0 

b  -  A  -  A  -  A  -  A  -  A  »  6AU 


or 

b  „  [6,  1,  0,  1,  1,  1,  1,  0] 

tne  values  of  the  basic  variables.  Another  way  to  compute  is,  of 
a  —  u  9 

course,  by  b  ■  b  -  0A  ,  the  usual  formula  for  updating  the  values 
of  the  basic  variables  is  the  simplex  method. 


-21- 


REFERENCE 


[l]  Dantzig,  George  B.  ,  and  R.  M  Van  Slyke,  "Generalized  Upper 

Bounding  Techniques  -  i,  "  ORC  64-17  (RR),  Operations  Research 
Center,  University  of  California,  also  to  appear  in  the 
Proceedings  of  the  IBM  Symposium  on  Combinatorial  Methods 
(March  1964). 


-22- 


